Not to appear in Nonlearned J., 45. 



A Very Close Binary Black Hole in a Giant Elliptical Galaxy 3C 

66B 

and its Black Hole Merger 

Satoru Iguchi 1 , Takeshi Okuda 2 , and Hiroshi Sudou 3 

ABSTRACT 

Recent observational results provide possible evidence that binary black holes 
(BBHs) exist in the center of giant galaxies and may merge to form a supermas- 
sive black hole in the process of their evolution. We first detected a periodic 
flux variation on a cycle of 93 ± 1 days from the 3-mm monitor observations of 
a giant elliptical galaxy 3C 66B for which an orbital motion with a period of 
1.05 ± 0.03 years had been already observed. The detected signal period being 
shorter than the orbital period can be explained by taking into consideration the 
Doppler-shifted modulation due to the orbital motion of a BBH. Assuming that 
the BBH has a circular orbit and that the jet axis is parallel to the binary angular 
momentum, our observational results demonstrate the presence of a very close 
BBH that has the binary orbit with an orbital period of 1.05 ± 0.03 years, an 
orbital radius of (3.9 ±1.0) x 10~ 3 pc, an orbital separation of (6.1 jlag) x 10~ 3 pc, 
the larger black hole mass of (1.2tol) x !0 9 M Q , and the smaller black hole mass 
of {7.0ttl) x 10 8 M & . The BBH decay time of (5.1^° 5 5 ) x 10 2 years provides 
evidence for the occurrence of black hole mergers. This Letter will demonstrate 
the interesting possibility of black hole collisions to form a supermassive black 
hole in the process of evolution, one of the most spectacular natural phenomena 
in the universe. 

Subject headings: black hole physics — galaxies: active — galaxies: elliptical and 
lenticular, cD — galaxies: formation — galaxies: individual(3C 66B) 
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Introduction 



IThe major galaxy mergers in galaxy formation inevitably lead s to the presence of a 



1980 



bin ary b lack hole (BBH) located in the center of a giant galaxy (IBegelman et al 
Yul 120021 ). The following observ ational results provide possible evidence of the presence 



of BBHs: double comp act cores ([Rodriguez et al.l 120061 ). two broa d-line emission system s 
( jBoroson fc Lauerl 120091 ). th e X-shaped morphology of radio lobes (jMerritt fc Ekerd 120 02). 



wiggl e patterns of radio jets (IRoos et al.lll993 



Abraham fc Carrara 



200 5]), periodic varia tions and periodic outbursts (jSillanpaa et al. 



1998 



1988 



Lobanov Sz Roland 



Valtaoia et al. 



2000 



Valtonen et al.ll2008l ). and the orbital motion of the compact core ( iSudou et al.ll2003l ). Strong 
evidence for BBHs as linking some observational results, however, has not yet been found. 
Clarifying the mechanism fo r the formation of BBHs is essential in the study of galaxy 
mergers in galaxy formation (ISanders et al.l Il988l ; iTreister et al.l 120101). in the role of black 
hole mergers in the evolution of supermassive black holes feaker et al. 2006), and in the 



Yu 


2002; 


Baker et al. 


2006) 



3C 66B is a near by low-luminosity F R I radio galaxy and giant elliptical galaxy with a 
redshift (z) of 0.0213 (IHuchra et al.lll999l ). and has the strong jet and counter jet that ex tend 
to abo ut 100 kiloparsecs (k pc) which have been observed at radio (IHardcastle et al.lll996l). in- 
frare d ( Tansley et al.|[2000 ). optical ( Macchetto et alll991 ) and X-ray waves ( Hardcastle et al 
200ll ). The orbital motion of the compact core in 3C 66B, which has a period of 1.05 ± 0.03 
years, was observed with a position accuracy of 10 mic ro arcseconds (//as ) by using phase- 
referencing Very Long Baseline Interferometry (VLBI) ( ISudou et al.ll2003l ). 



Several observational and theoretical studies have shown that the periodic flux variation 
and periodic outburst activities in the radio, optical, X-ray, and 7-ray light curves are asso- 
ciate d with the orbital mot i on of a BBH in the center of the galaxies ( iRieger fc Mannheim 
20001 ; |Pe Paoils et al.l 120021 . 12004 ). Here we report the detection of a signal periodicity in 
light curves from the compact core of 3C 66B, which indicates the presence of a very close 
BBH in the center of this object, and also present evidence for black hole mergers. We 
use a Hubble constant (H ) of 71 km s _1 Mpc -1 , the matter density (flu) of 0.27 and the 
vacuum energy of 0.73 in this Letter, resulting that an angular size or separation of 1 
milliarcsecond (mas) corresponds to 0.436 pc at the distance of 3C 66B. 
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2. Observations, Data Reduction, and Results 

The millimeter-wavelength flux variation for 3C 66B was monitored every two weeks 
from the end of 2003 through to the beginning of 2005 at 93.716 GHz using the Nobeyama 
Millimeter Array (NMA) of the National Astronomical Observatory of Japan (NAOJ), and 
every four weeks from the end of 2005 through to the middle of 2006 at 86.2 GHz using 
the Plateau de Bure Interferometer (PdBI) of the Institut de Radioastronomie Millimetrique 
(IRAM). These flux monitor observations using two- type millimeter arrays will enable us to 
eliminate any dependence on the characteristic features of each array, giving us more reliable 
and accurate results. 

In the NMA observation, 3C 84 and 0133+476 were observed as a visibility calibrator. 
Also, the flux densities of each calibrator were precisely derived by observing Uranus and 
Neptune. Phase fluctuations in the observation are caused by short-term variations of the 
excess path length in the atmosphere, which is mainly due to the time variation of water 
vapor pressure. Since these phase variations cause decoherence in the visibility data, it 
is necessary to correct this loss. The decoherence factors at all observation epochs were 



estimated from the synthesized maps of each visibility calibrator ( iKohno et al.l 120051 ). To 
accurately derive the flux densities of the observed object, the flux data at epochs with a 
decoherence of more than 20 % were flagged. 3C 454.3 and 0420—014 were observed for the 
bandpass calibration. The weighted mean of the flux density data that were estimated from 
each visibility calibrator was plotted. In the PdBI observation, the flux densities of 3C 66B 
were corrected by observing 3C 84 and the compact Hn region, MWC 349 and/or CRL 618, 
except 3C 84 on January 14 and MCW 349 on July 23 in 2006. The phase and bandpass 
calibrations were performed by 3C 66B itself. 

We first investigated the contribution of large-scale jet emissions into observed flux 
densities. A map of 3C 66B with millimeter jets can be made by combining the visibility 
data obtained from the NMA monitoring observations for all epochs (see Figure la). The 
total flux density of the millimeter map of 3C 66B at 93 GHz exceeds 500 mJy, while the 
sum of the expanding jet and counter jet is less than 24 mJy. We made the Spectral Energy 
Distributions (SEDs) of the jet and the counter jet from our data and other data at different 
wavelengths, respectively. We found that these SEDs follow (agree with) a simple power-law 
synchrotron model, even though the flux densities at all wavelengths were observed at the 
different dates (see Figure lb). This fact indicates that the observed flux variations of 3C 
66B were dominated by the unresolved core, not by the expanding jet and counter jet. 

Figure 2 shows the 3- mm peak flux monitor data obtained from the unresolved core of 3C 



(Lombl 


1976; 


Scarele 


1982) 



-4- 



results, we detected a periodic flux variation on a cycle of 93±1 days with a maximum-to- 
minimum amplitude ratio of 1.19 ±0.03, and non-periodic flux variation from the unresolved 
core. The variability timescale of the observed signal is shorter than that of the orbital 
period of about 1 year. Due to the orbital motion around the center of gravity in a BBH, the 
Dopp l er factor for the emitt i ng re gion is a periodical function of time ( IRieger fc Mannheim 



pp. 

2000l ; |Pe Paoils et al.l 120021 l2004j ). and the observed signal periodicity has a geometrical 



origin due to Doppler-shifted modulation. Comp ared to the intrinsic period, the ob served 
signal period is shortened by the relativistic effect (ICamenzind fc Krockenbergerlll992l ). It is 
also noted that the non-periodic flux variation may be caused by the non-thermal radiation 
in the flaring state that is emitted by a relativistic emission region not directly affected by 
a BBH (e.g., a shock in the knots of jet in the core that is not resolved with the angular 
resolutions of NMA and PdBI). 



3. Discussion and Summary 



Both the observed flux variation (see Figure 2) and the orbital motion (ISudou et al. 

20031 ) were well fitted by a simple model of sin function. This result indicates that the 



difference between the apogee and perigee velocities will be small, and that the Doppler 
boosting factor will be small. If they were large, the profile of the observed flux variation 
would be that of a spike like the periodic outburst activities in the flaring state. This 
indicates that the binary system in 3C 66B has a nearly circular orbit. As a consequence of 
the binary orbital revolution around th e center of gravity, the obse rved signal period (P bs) 
is related to the orbital period (Pk) by ( IRieger fc Mannheim! 120001 ) 



P obs = (1 + z) (1-/3 cos i)P k , (1) 

where (3 is v/c with the jet outflow velocity of v and the light speed of c, and i is the 
inclination angle between t he jet axis and the line of sig ht. The apparent speed in the 
relativistic jet is written by (IGinzburg &: Syrovatskiilll969h 



app 



(/3sinz)/(l — /3cosi) 



(2) 



Given the detected signa l period (P bs) of 93±1 days, the orbital period (pj of 1.05 ± 0.03 
years (ISudou et al.l 120031 ) . and the derived apparent speed (/3 ap p) of 0.30±0.04 (see Figure 
3), the following physical parameters for an inner-jet outflow that has the effect of the 
BBH orbital motion found in 3C 66B were estimated from equations ([T]) and (F2J); a speed 
(P=v/c) of 0.77±0.01 and an inclination angle (i) of 5.3±0.8 degrees in accordance with the 
standard model for jet bulk motion. Consequently, the derived Doppler factor (5) of 2.8 
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[=(1 — /3 2 ) 1 / 2 /(l — (3cosi)] is a small value that does not have a strong effect on the profile 
of the periodical flux variation. 

The jet is presumably linked to the accretion disk around a black hole. At least one of 
the two black holes in a binary system will have an accretion disk. If both of the black holes 
have an accretion disk, the relationship between the observed jet and the two black holes 
is not simple. Assuming that the observed jet is dominated by the accretion disk around 
either one of the two black holes, we show two schematic geometries of a BBH in Figure 4. 
The orbital motion velocity (uy) around the center of gravity can be written as u}Rsm(ut), 
where u = is the intrinsic orbital frequency. The observed flux modulation by Doppler 

boosting can be written by 

S{y) = 5(i) 3+Q S», (3) 

where S*{y) is the intrinsic spectral flux density in the comoving frame, a is the source 
spectral index (S oc and S(t) is the Doppler fa ctor including the time varia tion due to 

the orbital motion of the BBH, which is written by (IDe Paoils et al.ll2002l 120041 ) 



6{t) 



y/1 ~ (I 2 + /3 bl (t) 2 ) 

— (/3 cos i + /?bi(t) sin i) 



(4) 



where fiu(t) is Vb\(t)/c. From equation 
can be expressed as 



the observed maximum-to-minimum flux ratio 



3+a 



(5) 



By applying this model to the sin curve of the periodic flux variation caused by a BBH (see 
Figure 2), an orbital motion velocity (vu) of the BBH can be derived from equations (JH) and 
©, and an orbital radius (R) of (1.2 ± 0.3) x 10 16 cm [(3.9 ± 1.0) x 10~ 3 pc] can then be 
estimated when a—0. Although an upper limit for the orbital radius coul d only be estimated 
from the core motion using phase-referencing VLBI ( iSudou et al.l 120031 ). the orbital radius 
of a BBH can be directly derived by adding the periodic flux monitor data that we found 
successfully. 

From two types of geometry of a BBH system shown in Figure 4, an orbital separation 
(d) of two black holes is R(l + q) (see Figure 4a) or R(l+q)/q (see Figure 4b) when q = m/M, 
where m and M denote the mass of the smaller and larger black holes, respectively. Since 
the revolution period of a two-body problem with a circular orbit has a relational expression 
of Pk = 2irG~ 1 ^ 2 d 3 ^ 2 (M + m) -1 / 2 (where G is the gravitational constant), the mass of the 
larger black hole (M) is estimated to be 



{2u) 2 P^G- 1 R\l + qf 



(6) 
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or 



(27r) 2 P k - 2 G- 1 J R 3 (l + g) 2 /g 3 . 



(7) 



Ferrarese fc Merrittl (120001 ) demonstrated a tight co rrelation between the black hole masses 
and the velocity dispersions of their host bulges, and lMerritt fc Ferraresd (120011 ) then derived 
the relationship between central stellar velocity dispersions a nd black hole masses. Based 



on th is relationship, an estimated mass of (1.9±J;g) x 10 9 M flNoel-Storr. Baum. fc O'Deal 



20071 ) is adopted as the total mass of the two black holes [M(l + q)) for 3C 66B, where 
M Q is the mass of the Sun. From equations (jSJ) and ([7j), the following basic BBH physical 
parameters were estimated; the larger black hole mass of (1.2j:{j;i) x 109 M ©> the smaller 



black hole mass of (7.01^4) x 10 8 M , and the separation (d) between the black holes 



of (1.9 ± 0.3) x 10 16 cm [(6.I+09) x 10 3 pc] that corresponds to about 53+|; times the 
Schwarzschild radius (R s = 2GM/c 2 ) of (3.5±JJ) x 10 14 cm [(l.llju) x !0~ 4 pc, 23t£° AU]. 
This estimated separation is much short er than that of a celestia l object recently identified 
as being a potential sub-parsec BBHs by lBoroson fc Lauerl (120091 ). revealing the presence of 
a very close BBH in the center of 3C 66B. 

The decay time of a BBH in the center of 3C 66B due to the gravitational radiation is 



estimated to be (S.ll^ 5 ) x 10 2 years using the following equation (IPetera Il964j ; lYul 120021 ) 



5.8 x 10 



M 



V 10 S M P 



Q 



-3 



d 



O.Olpc/ q(l + q) 



The estimated decay time is much shorter than the Hubble time. The model of galaxy 
formation based on ACDM indicates that there were a large number of galaxy mergers, 
which led to giant elliptical galaxies. Centaurus A (CenA), which is the nearby radio galaxy 
and bright elliptical galaxy, has the pr ominent warped dust lane and gas disk in the central 



Baade fc Minkowski 



1954 ), suggestin g that Cen A experienced a recent merging 



region (e.g., 

event during a few 10 8 years (jStruve et al.ll2010l and referen ces therein). 3C66B is a giant 
elliptical g alaxy within the cluster Ab ell 347 ( IFanti et al.lll982l ) and has the prominent face-on 
dust lane (IVerdoes Kleiin et al.l 19991) . These indicate that 3C 66B differs from quasars with 
active merging ( iVolonteri et al.l 120091 ) and is one of nearby radio galaxies that experienced 
the galaxy merging event. Within a timescale of 10 8 years, two black holes caused after 
the galaxy merging event would come closer to each other due to the dynamical friction 
and gas interaction, a nd then merge into one due to the gravitational- wave radiation (e.g., 
Begelman et al.lll980l ). In conclusion, this shows a possibility that a very close BBH in the 
center of 3C 66B exists just before its merger into one black hole. The detection of strong 
gravitational-wave radiation is expected during the final stage of a black hole merger. 

In this Letter, the following assumptions are adopted; 1) the derived apparent speed 
of the jet is equal to the bulk flow speed; 2) the jet axis is parallel to the binary angular 
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momentum and is perpendicular to the accretion disk around a black hole by which the 
observed jet is formed; 3) the BBH has a binary orbit with an orbital eccentricity of zero 
(this is a circular orbit); and 4) the spectral index is zero at around our observing frequencies. 
It is theoretically suggested that the pattern speed is slower than the bulk flow speed of the 
jet (ILind fc Blandfordlll985l ). while the study on the relativistic bulk moti on in about 100 



AGN s demonstrates that there is no clear distinction between these speeds (iGhiseUini et al. 



1993!) . For instance, the ratio of these speeds for GRS1919+105 was estimated to be 0.8 
( jBodo fc Ghisellini 1995 ). As our estimated proper motion has an error of 15 % that is 
comparable to this uncertainty, the assumption 1) can be considered reliable in this Letter. 
The complexity of mergers and active galactic nuclei (AGNs) feeding may not support the 
assumptions 2) and 3). To evaluate these assumptions, it is necessary to monitor the flux 
variation with higher precision as closely as possible. If these assumptions are not entirely 
inappropriate, however, our derived physical parameters are not significantly changed. Since 
the synchrotron radiation from the AGN core is expected to be optically thin or peak at 
around our observing frequencies, we adopted the assumption 4) in this Letter. When 
a = 0.5 as compared with a = 0.0, the derived physical parameters of a BBH in 3C 66B 
are changed at a factor of 1 to 3. Using the Atacama Large Millimeter/sub millimeter Array 
(ALMA) including the Atacama Compact Array ( AC A) that is used t o obtain precise images 
of celestial objects with high angular resolution (jlguchi et al.l 120091 ). the mult i- wavelength 
profile of the periodic flux variation and periodic outburst activities of candidate objects 
can be monitored with precision of one mJy or better. The spike profile of the periodic flux 
variation will provide us with the orbital eccentricity of a BBH in 3C 66B and the relationship 
between the jet axis and the binary angular momentum. In addition, the multi-wavelength 
profile will improve our derived basic physical parameters of the BBH with high precision; 
the orbital radius, the orbital separation and the masses of the two black holes, and the 
decay time. 
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Fig. 1. — (a) The NMA map of 3C 66B at 93 GHz with the jet and counter jet with the 
synthesized beam of 3'.'0x2'.'8 with a P.A. -40 degrees, with contours at 1.1 x (-3, 3, 4, 5, 6, 7) 
mJy per beam, and with a peak flux of 126 mJy per beam, (b) Spectral Energy Distributions 
(SEDs) of the jet (filled circle) and the counter jet (open square) in the frequency range of 
radio to X-ray in 3C 66B. The flux densities at several wavelengths for the region of the jet 



and c ounter jet out to 8 arcsec from the core have been tabulated so far ( jHardcastle et al. 



20011 ). We estimated an integrated jet flux density of 21±1 mJy and an integrated counter- 
jet flux density of 2.5±1.8 mJy at 93.716 GHz, and additionally plotted these new points in 
the spectrums of the jet and that of the counter jet, respectively. 
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Fig. 2. — (a) Peak flux monitor data for the millimeter-wavelength core of 3C 66B with 
NMA (filled squares; 93.716 GHz) and PdBI (filled circles; 86.2 GHz), (b) Lomb-Scargle 
periodogram from these flux monitor data. We assume that the spectral index of the un- 
resolved core is zero between these two observing frequencies. The dotted horizontal lines 
present the false detection limit with the probability of 99 % and 95 %, respectively, on the 
presumption that the noise is white noise, showing the claimed periodicity of 93 ± 1 days 
with a 98 % probability. Shown in (a), the claimed periodic flux variation with a maximum- 
to-minimum amplitude ratio of 1.19 ± 0.03 was derived by using the minimum Chi-square 
model fitting between the observed data and the sin and cubic functions. 
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Fig. 3.— Four 8-GHz VLBI maps of the jet of 3C 66B observed using the Very Long 
Baseline Array (VLBA) of the Nati onal Radio Astron omical Observatory (NRAO) between 
2001 March 13 and 2002 June 14 (jSudou et al.ll2003l ). These maps were obtained using a 
restored synthesized beam of 0.8x0.8 mas with contours at 3a x (y2, 2, 2 a/2, ...) where a is 
0.5 mJy per beam, and were fitted with one Gaussian component in the core and three point 
components in the jet. The proper motions of the three jet components (dashed lines) were 
determined to be 0.21±0.03, 0.36±0.04 and 0.63±0.18 mas per year by using the five maps 
(including the map obtained on 2002 February 21; this map has been removed in this figure, 
however, because the date is close to that of the map on 2002 February 8). Considering a 
relativistic emission region directly affected by a BBH, we, in this Letter, use the proper 
motion of the most inner jet, giving an apparent speed (/3 ap p) of 0.30±0.04. 



13 




Fig. 4. — Two schematic geometries of a BBH in the elliptical galaxy 3C 66B. We assume 
that the BBH in 3C 66B has a circular orbit, that the jet is linked to the accretion disk 
around one of the two black holes, that the non-thermal radiation propagates outwards from 
the core along the jet with a velocity of v , and that the jet axis is parallel to the total angular 
momentum of the binary. The observed jet is formed by either (a) the smaller massive black 
hole (with a mass of m) or by (b) the larger massive black hole (with a mass of M). 



